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Abstract 

In the model of quintessino as dark matter particle, the dark matter and dark energy are unified 
in one superheld, where the dynamics of the Quintessence drives the Universe acceleration and its 
superpartner, quintessino, makes up the dark matter of the Universe. This scenario predicts the 
existence of long lived f as the next lightest supersymmetric particle. In this paper we study the 
possibility of detecting f produced by the high energy cosmic neutrinos interacting with the earth 
matter. By a detailed calculation we find that the event rate is one to several hundred per year at 
a detector with effective area of lkm 2 . The study in this paper can be also applied for models of 
gravitino or axino dark matter particles. 
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I. INTRODUCTION 



Recent astronomical observations strongly support a concordance model of cosmology, 
where the bulk of the content of the Universe is comprised of an unknown dark sector, of 
which about 23% is the pressureless dark matter and about 73% is dark energy with negative 
pressure that drives the present acceleration of the Universe. In the literature there are a 
lot of interesting and compelling models of particle physics which provide candidates for 
dark matter and dark energy. However, it is always much more desirable that a single 
model explains both for the dark sector. Recently, two of us (X.B. and X.Z.) with M. Li[l| 
have considered a class of Quintessence models, originally designed for the understanding 
of the current acceleration of the Universe and showed that after supersymmetrization the 
fermionic superpartner of the Quintessence, the quintessino, serves as a good candidate 
for the dark matter particle. Very much like the quarks and leptons in the same group 
representation in the grand unified theories, the components of the dark matter and the 
dark energy belong to one superfield in this model. In this scenario, the quintessino is 
the lightest supersymmetric particle (LSP) and the next lightest supersymmetric particle 
(NLSP) can be either neutralino or slepton, typically the stau. During the evolution of 
the universe, the NLSP freezes out of the thermal background, then decays into quintessino 
which makes up the present dark matter. 

Generally the quintessino dark matter particles interact with the ordinary matter very 
weakly, so the traditional techniques for the direct or indirect detection of the neutralino 
dark matter would not be useful for detecting quintessino dark matter. There are, however, 
several silent features of the quintessino dark matter model which make it distinguishable. 
Firstly, we should note that being one superfield the quintessino has the same interaction 
with matter as the quintessence and thus be severely constrained. As we know that the 
primary motivation to introduce the quintessence field is to drive the acceleration of the 
Universe, so we have to preserve the properties of quintessence after introducing interactions 
between quintessence and the matter. In order to keep the quintessence potential flat and 
avoid the long range force induced by quintessence we impose a global shift symmetry for 
its interactions Q], i.e., interactions which are invariant under Q — > Q + constant. One 
possible interaction is proposed in Ref. |2| £q 77 = jS-QF^F^ , which can be taken to 

pi 

test the quintessence model by observing the rotation of the polarization plane of the light 
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from distant sources. Another interesting form is £q// = j^-d^QJ 11 with J 11 the baryon or 

pi 



baryon minus lepton current [3] which gives rise to a new mechanism of the baryogenesis. 

Secondly, the BBN observation indicates that the light element 7 Li may be underabundant 
compared with the theoretical estimate j4|. Decaying particles after BBN might provide one 
way to solve the problem J] . As shown explicitly in Ref . (and argued for the gravitino 
dark matter models in Ref. P]) the electromagnetic energy associated with the non-thermal 
production of quintessino dark matter particles can play the role to reconcile the observations 
and the theory. 

Thirdly, the quintessino dark matter in this model is produced nonthermally. The prop- 
erty of the quintessino dark matter is characterized by the comoving free streaming scale. 
Depending on the time when the NLSP decay and the initial energy of the quintessino when 
it is produced, it can be either cold or warm dark matter. In the latter case it helps solve 
the problems of the cold dark matter on subgalactic scales p. 

Fourthly, this scenario for dark matter predicts the existence of long-lived NLSP with life 
time 10 5 — 10 8 sec. For the stau NLSP it can be produced and collected on colliders|Q] and 
the properties of quintessino can be studied by examiningthe stau decay. The scenario with 
gravitino LSP and stau NLSP is studied in the literature [7]. For the case of quintessino LSP 
the stau will have different decay modes. 

In the present paper we are not attempting to study all the features of quintessino dark 
matter mentioned above, instead we will focus on the study of the NLSP stau production 
and detection. We leave other phenomenological studies for future publications 

The cosmological stau will not survive until today, however it can be produced in the ex- 
tremely high energy astrophysical processes. Experimentally these stau can then be observed 
by the large cosmic particle detectors on the earth, such as the L3C|9], SuperKamiokande . 
or several proposed neutrino telescopes^]. Since the stau are much heavier than other sta- 
ble or long lived charged particles, such as the electron or the muon, they are expected to 
leave distinct trajectory in the detectors. 

There are different sources of the stau fluxes. However, since the stau has life time about 
one year we expect the stau flux from distant objects, such as AGNs, should be quite low. 
Therefore, most stau signal will come from the earth. One stau flux comes from the high 
energy cosmic neutrino collision with the earth matter and the other one comes from the 
collision of the high energy cosmic proton and nuclei with the atmosphere. In this paper 
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we will present a detailed calculation of the stau flux produced in the simpler case by the 
cosmic neutrino collision with the earth matter and leave the discussions on the atmospheric 
flux of stau in another workQ]. In a recent paper by Albuquerque, Burdman and Chacko 
the production of the NLSP stau by the high energy neutrino flux has been proposed as a 
possibility of probing for the supersymmetry breaking scale using the neutrino telescope[12]. 

The stau flux thus generated depends on three ingredients: the flux of the primary cosmic 
ray of the incident neutrinos, the cross sections of the neutrino and nucleon and the range 
of the charged stau on the earth. In this paper we will show that there will be one to several 
hundred of stau produced per year at a detector with effective area of one km 2 . The paper 
is organized as follows: in Sec. II, we study how quintessino dark matter is produced by the 
stau NLSP decay. In Sec. Ill, we calculate the neutrino-nucleon interacting cross sections. 
In Sec. IV we discuss the energy distribution of f and its range. In Sec. V, we describe our 
calculations and present the results. Finally, Sec. VI is our conclusion. 



II. QUINTESSINO PRODUCTION FROM NLSP STAU DECAY 

In Ref. Q| we focus our study of quintessino dark matter on the case of neutralino being 
NLSP. Therefore we will first study the viability of stau NLSP to produce the quintessino 
dark matter in this section. 

The most general form of the interaction between the quintessence and the leptons that 
obeys the shift symmetry is given byQ] 

A?// = )*M<§fiitffiR + <'■,(, ,-".(;,,:) , (i) 

with A being the cutoff scale and c^- being coupling constant in the effective theory. The 
astrophysical and laboratory experiments put a lower bound on A of about 10 w GeV [lj]. 
Here we take a simple form of the interaction between quintessence and the lepton r 

C-QSi = £ WtVt ■ (2) 

After supersymmetrization and using the on-shell condition of r we get the relevant coupling 
for our discussion 

C Qrf = ^-L(QP LT r R + QP R rfl) + h.c. , (3) 
where f L and tr are the left- and right-handed stau respectively. 



The quintessino is produced by the stau decay, r — * Q + r. From Eq. (j3J), the time scale 
for the stau decay is given by 

. ~ 10 6 sec / A \ 2 flTeV\ n 2 , 2 

t{t^Qt )«——(- 1-4 2 , (4) 

<r \A GUT / Km?/ 

where xq = m^/m- and Aqut ~2x 10 16 CreV being the scale for grand unification. Taking 
10 10 CreV ^ A < 10 13 CreV, the stau decays before BBN. In this case there will be no 
constraints on the process from BBN. 

For higher scale A, the stau decays after BBN. The electromagnetic energy released in the 
stau decay may de Stroy the successful predion ou the h g ht ahuudauce of BBN0 

and the Planckian spectrum of the cosmic microwave background (CMB)Q). Therefore we 
have to consider the constraint from BBN and CMB. The electromagnetic energy released 
in the stau decay can be written as 

1 - x 2 

tEM = e EM N f = 0.75 x 10- g GeV & , (5) 

where eem is the initial electromagnetic energy released in each stau decay and N f = 
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TeV 






m Q . 





-Q2£ I is the number density of stau normalized to the number 



density of the background photons 0. As the argument given in Ref. 0, we take €em = \E T 
for simplicity, where E T is the energy of r in the decay products. 

In Fig. Q we plot the decay time and energy release for m f from "iOOGeV to 3TeV and 
xq from 0.1 to 0.9. The shaded region is excluded by BBN. The best fit point (t,£em) — 
(3 x 10 6 sec, 10~ 9 Ge\^) is covered in the parameter space, which can suppress the level of 7 Li 
and make it to be consistent with the observation 13j. 

We will not further consider the constraint from BBN on the stau three-body hadronic 
decay in this work. Compared with the recent analysis given by Feng et. azjl3|, we find 
there is viable parameter space around the best fit point even after imposing the hadronic 
constraint in our model. Therefore, we get the conclusion that the quintessino dark matter 
produced by stau decay is a viable dark matter model, which is consistent with the present 
experimental constraints. 
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FIG. 1: Lifetime r and energy release £em i n stau decay for rrif from 300Gey to 3TeV (from 
right to left) and xq from 0.1 to 0.9 (from top to bottom). We take c = 1 and A = 2 x 10 16 GeV. 
The contour of the chemical potential for photon distribution function is \i = 9 x 10 -5 , with the 
region to the right of it being excluded by CMB 15]. The shaded regions are excluded by BBN 
data 13]. The circle represents the best fit region with (t,£em) = (3 X 10 6 sec, 10~ 9 GeV) 13]. 



III. SUPERSYMMETRIC NEUTRINO-NUCLEON CROSS SECTION 

Probing for the long lived stau provides an indirect support for our dark matter model 1 . 
In this section we calculate the inclusive cross section of stau production from the collision 
of high energy cosmic neutrinos on the earth matter. 



1 In the case of the gravitinoQ or the axino J8j as the LSP and the candidates for dark matter particle, 
stau can serve as the NLSP in the similar way. Our studies in this paper can be easily generalized into 
these cases. 
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FIG. 2: The cross sections of the inclusive process Eq. © for Mi = Mi = fj, = 250GeV, = 
m eR = 150GeV, and uig = m u ~ R = m d ~ R = 150, 300, 500, 800Gey . The curve with the smallest cross 
section is for M\ = Mi = [i = AOOGeV, = m e ~ R = 250GeV, and mQ = m u ~ R = m d ~ R = 800GeV. 
The cross section for the standard model charge current process is also plotted in the figure. 

The process which we are interested in is 

u + N^l + q + X^ 2f + X' , (6) 

where N = |(n + p) is an isoscalar nucleon, X (X') refers to any type of hadrons and 
leptons. At the parton level the corresponding process is v + q — > I + q with the exchanges 
of the chargino, x + 1 01 neutralino, x° 111 t ne i-channel, where q stands for the valence and 
sea quarks. The I and q in Eq. (JHJ) will quickly decay into two NLSPs. For the calculation 
of the inclusive cross section in Q, we use the automatic Feynman Diagram Calculation 
, which adopts the CTEQ6-DIS parton distribution function [20J for the 

isoscalar nucleon. 

For the soft SUSY breaking terms we take the parameters as = mg fl = 150GeV in the 
calculations, which gives rise to a lighter stau of about 143GeV. The gaugino masses and 
Higgsino mixing parameter are taken as Mi = M2 = /i = 250^6^, so that the neutralinos 
and charginos are heavier than f . These parameters give the lightest neutralino mass m x o = 
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181GeV and the lighter chargino mass m x + = 193GeV. In Fig. |21 we plot the cross sections 
for mQ L = rriq R = 150, 300, 500, 800GeV, which are the most sensitive supersymmetric 
parameters to the cross section. Another set of the soft parameters are taken as = 
m iR = 250GeV, M x = M 2 = /i = AOOGeV and m§ h = rriq R = 800GeV. (We will denote this 
set of parameters by mj = 250GeV in the text.) In the figure, the inclusive cross section 
of the charged current process + iV — > //+ anything is also plotted. One can see that 
the supersymmetric cross section is about 3 orders of magnitude smaller than the standard 
model one. 



IV. THE DISTRIBUTION OF NLSP ENERGY AND ITS RANGE 

When the staus propagate in the earth they lose energy through ionization and radia- 
tion processes: bremsstrahlung, pair production, and photonuclear interactions. While the 
energy loss due to ionization is only slowly logarithmically increasing with the energy, the 
radiation processes cause a loss which in the high energy limit is proportional to the stau 
energy. The energy loss is usually expressed as 

dE 9 



AX " + l3Ef - (7) 

where a is due to ionization and ft is due to radiation. We take a = 2 x 10~ 3 GeV g~ 1 cm 2 , 
which is the same as the ionization loss coefficient of fi in the rock^sj. Since ft is proportional 
to the inverse square of the incident particle mass, we take ft = ft^ ■ (jji^) ■ For /i, ft^ = 
3.9 x 10~ 6 Q~ 1 cm 2 jl^jj. We then obtain the range of the stau with initial energy Eq 

fi(£„) = ilog^±|^, (8) 

ft a + ftrrir 

where we have assumed that the energy threshold of the detector can be very small. 

The range of the stau is actually a distribution as a function of its energy. For simplicity, 
we use the range for the average energy of stau, R((E?)), in the calculation of the stau flux. 
Since the collision is through t-channel, the stau from the slepton decay and that from the 
squark decay have very different energies. To discriminate them, we denote them by fx and 
f 2 , and the corresponding rang R fl and Rf 2 respectively. 

For fx, we assume it has the same energy as the initial slepton, since all the three flavor 
sleptons have the similar masses. While for f 2 , we have to calculate the energy distribution 
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FIG. 3: Angular (cos 9) and energy distribution of the final particles in the v + N CM system, 
with the beam direction referring to the direction of the incident nucleon. The figures represent 
the angular distribution of f\ and f 2 with the beam direction and the energy distributions for f±, 
T2, I' respectively and the angle between the two staus. We take the CM energy of v and N at 
1000 GeV. At higher energies the angular distribution of fi 2 is more concentrated at —1 and the 
angle between them becomes smaller, while the energy distribution is almost not changed. 

for the process q — > f% + V + q' . In our calculation the energy distributions for the process 
v + N — > f\ + f 2 + q' + V + X are obtained by using the Monte Carlo method to calculate 
the four-body final state. In Fig. 01 we plot the angular and the energy distribution of 
the final state particles of the process above in the v + iV CM system. The angle between 
the two stau is plotted in the last panel. From the figure we see that the two stau are 
almost in the same direction as the incident neutrinos. This is easy to be understood since 
the struck quark carries only small fraction of the total momentum of the nucleon. As 
the energy increases, the f\ and f 2 are more concentrated in the forward direction of the 
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incident neutrinos. Therefore, the transverse momentum of the slepton in the CM system is 
generally much smaller than the momentum of v and N. Considering that the momentum of 
the struck quark is widely distributed from just above the threshold to the total momentum 
carried by the nucleon, we find it is difficult to estimate the transverse momentum of the 
slepton. 

In the v + N CM system, the average energy of f\ changes from ~ 65%E* to ~ 75%E*, 
while the average energy of T2 changes from ~ 15%E* to ~ 35%E*, both depending not 
sensitively on the energy of the neutrino, E*, and the squark mass. We fix E~ = 70%E*, 
E~ = 20%E* in the CM systerm as a simplification of the calculation. Assuming that the 
two stau are both in the forward direction, we boost them to the laboratory system. We 
then get 

E fL ~ — ( E *f + \Pf\) « « 70%(20%)£„ , (9) 

1(2) m N m N 

where E u is the energy of the primary cosmic neutrinos, and we ignore the f mass in the 
second step, which is a rough approximation at low energies. 



V. EVENT RATE 
A. Neutrino flux 



Waxman and Bahcall (WB) have set an upper bound on the high energy diffuse neutrino 
flux from the observed cosmic ray flux at high energies. Depending on the evolution of 
the source activity, they gave the limit on muon neutrino and anti-neutrino extra-galactic 
fluxQ 

f^p] = (1 „ 4) x 10- 8 • ( l -^-\ 2 GeV^cm'h-hr- 1 . (10) 

\ dE /WB \ E J 

If "unknown source" of protons is taken into account, the low energy neutrino flux can 



be raised as dcj) v /dE oc E 2 [1 + 0.1((10 8 GeV) / E)} [2l|. However, the neutrino flux at the 



energies below about ~ lO^GeV is bounded by AMANDA experiment |22j|. i.e., d4> u /dE < 
8.4 x 10 -7 E~ 2 GeV cm" 2 s' 1 sr^ 1 . We refer to the maximal value of flux (fTTIjl as 'WB1' flux, 
while the flux raised at low energies and bounded by AMANDA experiment below ~ 10 6 GeV 
as the WB2' flux. 

The WB bound holds for the optically thin neutrino sources. Mannheim, Protheroe and 



Rache (MPR) extend the calculation to "optically thick" sources to the nucleons[23|. They 

10 



o 

FIG. 4: The path of neutrino flux penetrates the earth to the detector at the angle 9. 

give an upper bound of diffuse neutrino flux which is almost 30 times larger than the WB1 
limit in Eq. (ITUJ) . 

When the extra-galactic neutrinos propagate to the earth, we expect they contain the 
three flavor neutrinos in ratio of 1:1:1 due to the neutrino oscillation. The supersymmetric 
particle production rate should be independent of the neutrino flavor since we take all the 
sleptons almost degenerate. Therefore the three flavors of neutrinos should contribute with 
the same signal rate, except that v e has more strongly scattering effect with the earth 
electron. Here, we calculate the supersymmetric particle signals at the detector contributed 
by the + flux with WB1, WB2 and MPR limit. 

B. Calculation 

The earth becomes opaque for neutrinos above about AOTeV from the nadir due to the 
charged- current interaction between neutrino and nucleon. Therefore the scattering of z/ M to 
H has to be taken into account. 

In Fig. H]we show the path of a neutrino flux penetrating the earth before it reaches the 
detector located at the point O. We assume that the primary flux at the point A is <p UQ , 
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then the flux at the point x is given by 

Mx) = <P» e- J > (r) ^ dR , (11) 

where the integration is along the path from the point A to x, with n(r) the nucleon num- 
ber density, the standard model charged-current cross section and R is between and 



2R* cos 9, with = 6371 km being the radius of the earth. The number density n 



r) is 



given by n(r) = p(r)NA, where p(r) is the density profile of the earth given in Ref. 24] and 
Na is the Avogadro's number. 

Assuming the stau produced below the point B can reach the detector, we then get the 
flux of stau at the detector O 



MB)e-f y B n{r) ^ dR a f n{r)dy 
V {B)^ (l - e -JWir) f (12) 



where both integration are along the path and o> is the supersymmetric cross section. It 
should be noted that at the point corresponding to Rf x there is only one stau finally arriving 
at the point O, while at the point corresponding to Rf 2 , there can be two stau arriving at 
the point O. 

At large angles, the range of the stau may be larger than the actual depth of the earth 
that a neutrino has penetrated. In this case we take the actual depth the neutrino penetrated 
to calculate the probability that a stau can be produced by a neutrino. 

In Fig. 03 we plot the flux E^dfy/ ' dE v as a function of the incident neutrino energy 
for WB1, WB2 and MPR neutrino fluxes and for m~ q = 150, 300, 500, 800GeV. For m~ q = 
150Ge^ we show the results of WB1, WB2 and MPR fluxes. For other parameters we plot 
only the result of WB1 in the figure for clarity. We notice that the curves become almost 
flat at the energy above ~ 10 s GeV. The reason is that for the energies above ~ 10 8 CreV 
the range of f is so large that it is comparable with the radius of the earth. Therefore in 
most directions the <pf is proportional to the primary neutrino flux according to Eq. (|12j) . 
i.e., <p u (B) = <p uo for R? > R{6). 

Note that the contribution to the total f events comes mainly from the low energy neu- 
trino flux. (That will be very clear in a dN-J dE v — E v figure.) The reason is that the neutrino 
flux decreases with energy with the power law index of —2, while the supersymmetric cross 
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FIG. 5: Flux E^dfy/dE 1/ as a function of the incident neutrino energy for WB1, WB2 and MPR 
neutrino fluxes and rriq = 150, 300, 500, 800GeV . The flux for the set of soft parameters with 
mi = 250GeV is also plotted. Similar quantity for [i flux is plotted in figure. The three curves for 
[i flux calculated for the MPR, WB2 and WB1 neutrino fluxes respectively are shown from the top 
to the bottom. 







{mq =) 150GeF 


300GeF 


500GeV 


SOOGeV 


m l = 250GeV 


WB1 


1.8 x 10~ 4 


6.6 x 10" 5 


1.5 x 10~ 5 


3.3 x 10~ 6 


6.5 x 10~ 7 


5.8 x 10~ 7 


WB2 


3.2 x 10~ 3 


3.7 x 10~ 4 


4.9 x 10~ 5 


7.0 x 10~ 6 


1.0 x 10~ 6 


9.4 x 10" 7 


MPR 


5.4 x 10~ 3 


2. x 10~ 3 


4.4 x 10 -4 


1.0 x 10~ 4 


1.9 x 10~ 5 


1.7 x 10" 5 



TABLE I: Event rate of stau at the detector for per year per m 2 . 



section increases with energy (well above the threshold) with a smaller power law index of 
only about 0.4. 

The corresponding curves for the /i flux is also plotted in the figure for the WB1, WB2, 
and MPR neutrino fluxes. It is quite noticeable that the stau flux is even comparable with 
the /i flux in the case of irtq = 150GeV. 
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Finally, the total event rate for one year is given by 

event rate = lyear ■ 2, J frdE v d cos 0. (13) 

In table [fl we give the event rates per year per m for the WB1, WB2 and MPR fluxes for 
different squark masses. The total \x events per year per m 2 is also given in the table. It is 
quite interesting to notice that even in the most conservative case we may observe one such 
a event at a detector with lkm 2 effective area per year. 



VI. CONCLUSION 



In this paper we have considered the phenomenology of the model with quintessino as the 
lightest supersymmetric particle and the dark matter candidate [1]. In this scenario there 
exists a long-lived charged heavy particle, usually the lighter f. It is quite possible to detect 
such particles in cosmic ray detectors with large effective area. Discovery of the supersym- 
metric particles of cosmic source will be a valuable complementary for supersymmetry search 
at high energy colliders. 
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